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Abstract 

Data from the recent UH-60A Airloads Test in the National Full-Scale Aerodynamics 
Complex 40- by 80- Foot Wind Tunnel at NASA Ames Research Center are presented and 
compared to predictions computed by a loosely coupled Computational Fluid Dynamics 
(CFD) /Comprehensive analysis. Primary calculations model the rotor in free-air, but initial 
calculations are presented including a model of the tunnel test section. The conditions stud- 
ied include a speed sweep at constant lift up to an advance ratio of 0.4 and a thrust sweep 
at constant speed into deep stall. Predictions show reasonable agreement with measurement 
for integrated performance indicators such as power and propulsive but occasionally deviate 
significantly. Detailed analysis of sectional airloads reveals good correlation in overall trends 
for normal force and pitching moment but pitching moment mean often differs. Chord force is 
frequently plagued by mean shifts and an overprediction of drag on the advancing side. Loca- 
tions of significant aerodynamic phenomena are predicted accurately although the magnitude 
of individual events is often missed. 
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Figure 1 - UH-60A Airloads Rotor in 40- by 80- 
Foot Wind Tunnel. 


at matching conditions from the previous full-scale 
flight test and small-scale DNW wind tunnel test to 
assess rotor and wind tunnel scaling issues. Finally, 
unique slowed-rotor simulations were performed 
at reduced RPM (40% and 65%), achieving advance 
ratios up to 1.0. A complete enumeration of data 
acquired during the test is available in Ref. 1. 

Unfortunately, not all test data is perfect and no 
simulation is exact. The goal of this work is to pro- 
vide an initial correlation between measured data 
and a state of the art simulation. This correlation 
is intended to help discover flaws in experimental 
technique while at the same time identifying op- 
portunities to enhance rotorcraft simulation tech- 
nology. 


Introduction 

Testing was successfully completed in May 
2010 on a full-scale UH-60A rotor system in the 
USAF's National Full-Scale Aerodynamics Com- 
plex (NFAC) 40- by 8o-Foot Wind Tunnel [1]. 
The primary objective of this NASA/Army test 
program was to acquire a comprehensive set of 
validation-quality measurements on a full-scale 
pressure-instrumented rotor system at conditions 
that challenge the most sophisticated modeling 
and simulation tools. A secondary objective was 
to meet one of the original goals of the UH-60A 
Airloads program: to provide data to evaluate the 
similarity, or lack thereof, of measurements be- 
tween small-scale wind tunnel [2], full-scale wind 
tunnel (current test), and full-scale flight test [3]. 

The test hardware included the same rotor 
blades used during the flight test. Figure 1 shows 
these blades installed on the NFAC Large Rotor 
Test Apparatus (LRTA) in the wind tunnel test 
section. Key measurements included rotor per- 
formance, blade loads, blade pressures, blade dis- 
placements, and rotor wake measurements using 
large-field Particle Image Velocimetry (PIV) and 
Retro-reflective Background Oriented Schlieren 
(RBOS). 

Data were acquired over a wide range of test 
conditions, including speed sweeps at l-g sim- 
ulated flight conditions and parametric thrust 
sweeps (up to and including stall) at various com- 
binations of shaft angles and forward speed. These 
conditions included airspeeds up to 175 kt and 
thrusts up to 32,000 lb. Data were also acquired 


Methodology 

The complexity of rotorcraft aeromechanics is 
not easily modeled. The analytical results pre- 
sented herein were obtained using two separate 
codes — each a specialist in a particular aspect of 
rotorcraft simulation — loosely joined for this pur- 
pose. This section will describe the two codes and 
how they work together to maximize simulation 
accuracy and efficiency. 

CAMRAD II 

CAMRAD II belongs to a family of software 
known as "Comprehensive Codes" for the analysis 
of rotorcraft. These analyses incorporate a myriad 
of models to simulate the different aeromechani- 
cal subsystems of rotorcraft. CAMRAD II brings 
together a multibody dynamics model, a nonlin- 
ear finite elements structural model, and an aero- 
dynamics model based on lifting line theory [4]. 
CAMRAD II has seen a great deal of use in the 
simulation of the UH-60 aircraft in a variety of 
flight conditions [5-9]. 

The CAMRAD II structural dynamics model for 
the UH-60A has been decades in development by 
NASA and the U.S. Army. The specific model used 
here was refined by Yeo et al. in 2004 [10]. The 
model simulated the rotor using 7 l-D structural 
beam elements and 20 aerodynamic panels. The ro- 
tor was trimmed using Newton-Raphson iteration 
on collective and cyclic to meet specified trim tar- 
gets. 
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OVERFLOW 2 

All Navier-Stokes CFD analysis presented herein 
was performed using OVERFLOW 2 version 2.2b 
[11]. OVERFLOW 2 is an overset, structured-mesh 
flow solver developed at NASA. For two decades 
the OVERFLOW solver has served to analyze a 
variety of rotorcraft under a wide range of flight 
conditions [12]. OVERFLOW 2 offers a wide vari- 
ety of numerical schemes, turbulence models, and 
boundary conditions. For the present study, OVER- 
FLOW 2 was run with 4 th order central differencing 
and 4 th order artificial dissipation in space. Time 
marching was performed using a 2 nd order dual 
timestepping scheme. Turbulence was modeled 
near blade surfaces using the Spalart-Almaras one- 
equation model with rotational corrections. The 
turbulence model was deactivated in regions one 
chord length or further from the rotor blades to 
reduce numerical dissipation of the wake. Blade 
surfaces were modeled as viscous, adiabatic walls; 
outer boundaries were modeled using a character- 
istic condition imposing freestream quantities. 

OVERFLOW 2 computes the flowfield by dis- 
cretizing the Navier-Stokes equations on a series of 
overset, structured grids. Grids modeling the ro- 
tor blades were body-fitted and curvilinear. These 
grids, often called near-body grids, extended ap- 
proximately one chord length from the blade sur- 
face. The near-body grids were nested within one 
or more grids, called off-body grids, which filled 
the space between the rotor and the boundary of 
the computational domain. The OVERFLOW 2 
model included a notional hub, but the LRTA and 
wind tunnel struts were not modeled. All grids 
exchanged flow information in regions of overlap. 
The amount of this overlap was sufficient to sup- 
port full 4 th order accuracy at the boundaries. 

The grid surface is the latest definition derived 
from the as-built CAD model. Notable differences 
between this surface definition and that of prior in- 
vestigations [14,15] are an outboard shift of the trim 
tab by approximately 4 inches, a small bump on the 
upper surface near the blade grip, a slight thinning 
of the airfoil very near the tip, and a blunt trailing 
edge across the whole span. 

The near-body grid representing the bulk of 
each rotor blade had an "O" topology with 157 
points wrapping around the blade chordwise, 145 
points along the blade span, and 75 points nor- 
mal to the surface. The initial spacing at the blade 
surface had a y + value of 1. This grid system is 


illustrated in Fig. 2. The baseline near-body grid 
system (the hub and all four blades) contained 
approximately 10.7 million points. 

Two different off-body grid systems provided 
two different simulation scenarios: 

The first system used a series of ever larger 
Cartesian grids to create shells expanding outward 
from the near-body grid set. The grid point spacing 
within each shell is twice that of the shell imme- 
diately preceding it. The finest off-body grid had 
a spacing equivalent to 0.1 tip-chord lengths in 
all three directions. Seven such shells created a 
cubic computational domain spanning ten rotor 
radii in every direction. Domain boundaries in 
this scenario were set to a freestream characteristic 
condition thereby simulating a rotor operating in 
free air. Flow in the off-body grids was treated as 
inviscid and the turbulence model was deactivated. 
This off -body grid set consisted of 15.5 million grid 
points and is depicted in Fig. 3. 

The second off-body grid system used a sin- 
gle grid to envelop the near-body grid-set forming 
a computational domain that mimicked the size 
and shape of the wind tunnel test section. This 
grid made no attempt to model the contraction or 
expansion sections of the tunnel but rather main- 
tained the test section profile for 4.6 rotor radii 
upstream and downstream. The rotor plane was 
located 20.4 feet above the simulated tunnel floor. 
The inlet end of the tunnel grid employed the 
same freestream characteristic boundary condition 
as the free air system, and the outlet plane used 
a boundary condition that ensured conservation 
of mass within the tunnel grid. Viscous terms of 
the Navier-Stokes equations were discarded for 
this grid. This paradigm was initially proposed by 
Chang et al. [16] This off-body grid contained 4.2 
million grid points and is pictured in Fig. 4. 

In addition, a refined version of the free-air grid 
system was prepared to ascertain grid convergence 
of the simulation. This grid system reduced grid 
point spacing in both the near-body and off-body, 
approximately doubling the number of grid points 
to 50.7 million. 

Unless otherwise indicated, all predicted data 
presented herein were computed using the baseline 
near-body grid and the baseline free-air off-body 
grid set. This near-body and off-body combination 
totals 26.2 million grid points. To reduce computa- 
tion time CFD simulations were run using a hybrid 
distributed/shared memory scheme with 160 Mes- 
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Figure 2 - OVERFLOW 2 Near-Body Grid System. (Not every point shown.) 



Figure 4 - OVERFLOW 2 Wind Tunnel Off-Body Grid System. (Not every point shown.) 
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sage Passing Interface (distributed memory) ranks 
and eight OpenMP threads (shared memory) per 
MPI rank for a total of 1280 parallel tasks. OVER- 
FLOW 2 required approximately 31 minutes to ad- 
vance the solution for this configuration 14 of a ro- 
tor revolution using 1280 CPUs of an SGI Altix ICE 
computer. 

Coupling Methodology 

CAMRAD II uses a lower-fidelity aerodynamics 
model than that available in modern CFD codes, 
and most CFD codes lack the sophisticated Com- 
putational Structural Dynamics (CSD) and trim ca- 
pabilities of comprehensive codes like CAMRAD II. 
Coupling a CFD code (e.g. OVERFLOW 2) to a 
comprehensive code (e.g. CAMRAD II) marries the 
strengths of the two approaches and produces the 
highest-fidelity solution currently possible. For this 
study, coupling is achieved by alternate execution 
of OVERFLOW 2 and CAMRAD II. At the end of 
each code's turn to execute, it passes data to the 
next code. The data passed from OVERFLOW 2 
to CAMRAD II is airload data integrated from its 
Navier-Stokes model of the UH-60 rotor. This air- 
load data is used to augment CAMRAD II's inter- 
nal aerodynamics model (which consists of airfoil 
tables and a lower-order wake model). At the end 
of its execution, CAMRAD II generates updated 
control positions and a description of how the blade 
deforms elastically as it revolves around the shaft. 
These quantities are used to give OVERFLOW 2's 
grids a realistic motion in response to the aero- 
dynamic environment. This algorithm, called the 
delta coupling technique, was pioneered by Tung et 
al. [13] and implemented in OVERFLOW by Pots- 
dam et al. [14] Significantly improved airloads pre- 
diction capability has been demonstrated for the 
UH-60A rotor in steady level flight conditions us- 
ing this loosely coupled approach [14,15]. 

The CFD solution is advanced 14 revolution dur- 
ing a coupling iteration because this allows each 
of the rotor's four blades to sweep through a full 
quadrant of the rotor disk. Taken in aggregate, the 
four blades thereby determine the airloads at every 
azimuth for every coupling iteration. 

Convergence of the coupling process was deter- 
mined by monitoring blade airloads for periodic- 
ity. When the airloads did not vary significantly 
from one coupling iteration to the next, the solu- 
tion was judged to be converged. For the present 


analysis, this generally occurred after 24 coupling 
iterations. Since OVERFLOW 2 was allowed to it- 
erate for 14 revolution between coupling exchanges, 
this equates to 6 full revolutions for the converged 
solution. A fully converged coupled solution re- 
quired approximately 17 hours to compute for the 
baseline grid on 1280 SGI Altix ICE processors. 

Selected Test Conditions And Trim Ap- 
proach 

For the present investigation, two parametric 
sweeps were selected for analysis. First, a speed 
sweep was selected to test simulation accuracy over 
a wide range of advance ratios. Second, a thrust 
sweep was selected to test simulation accuracy un- 
der conditions ranging from a lightly loaded rotor 
through deep stall. 

The selected speed sweep ranged from /.< = o. 1 5 
to 0.4 with a constant Q\ J <7=0.09 and a constant tip 
Mach number of 0.65. During testing, lift, propul- 
sive force, and hub moments were trimmed to 
match nominal values at each flight condition by 
varying a s along with collective and cyclic pitch. 
In CFD calculations it is difficult to change the 
shaft angle once it has been set, which compli- 
cates its use as a trim control in adjusting lift and 
propulsive force. Instead, these flight conditions 
were matched in simulation by setting a constant 
a s as indicated by tunnel data while trimming Cj 
and hub moments to match tunnel values using 
collective and cyclic pitch. 

The selected thrust sweep was conducted at 
}i = 0.3 with M t i p = 0.625, a s = o°, and the rotor 
trimmed for minimum hub moments. The thrust 
was varied by changes in collective pitch up to 
stall. Absolute collective angles for a given thrust 
are generally different between simulation and test. 
To nullify this offset, a baseline case was selected 
at Cx/u=o.o8 and the simulation was trimmed to 
match the measured thrust. Collective deltas were 
derived from test data relative to this baseline and 
then applied to the simulated baseline to produce 
the remaining target points. 

In all cases, free-air simulations used corrected 
shaft angles, a c , derived from the geometric shaft 
angle, a s , by applying a Prandtl-Glauert wall cor- 
rection [1]. Simulations including tunnel walls used 
uncorrected shaft angle data. 
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Results 


The 2010 UH-60A Airloads Test is a valuable 
resource not only for the breadth of flight condi- 
tions obtained but for the depth of detail available 
in each measurement. Failure of analysis to pre- 
dict bulk loads accurately can be diagnosed by ex- 
ploring detailed measurements and how they com- 
pare to quantities extracted from simulation results. 
This approach will now be applied to the present 
simulation methodology first for the selected speed 
sweep and then for the selected thrust sweep. 

Speed Sweep 

Integrated Results 

Figure 5 plots the value of the trimmed thrust 
coefficient versus advance ratio for the range of 
speeds simulated. Since thrust was a trim target, 
it was expected that the computed points would 
lie very close to the measured data. Indeed this 
is the case for the data derived from CAMRAD II 
output. However, the data computed by integrat- 
ing the pressure and viscous forces acting on the 
CFD grid yield a somewhat higher thrust. The 
average difference between the thrust reported by 
CAMRAD II and OVERFLOW 2 is 2.5%. There are 
a number of factors contributing to this difference. 
First, small discrepancies in the local twist between 
the OVERFLOW 2 model and the CAMRAD II 
model can alter the computed thrust by altering 
in-plane and out-of-plane contributions by nor- 
mal and chord force. Second, the airloads passed 
from OVERFLOW 2 to CAMRAD II are only the 
loads integrated from the main blade grid and ne- 
glect any forces contributed by the tip and root 
caps. The coupling mechanism makes no effort to 
fully conserve forces between OVERFLOW 2 and 
CAMRAD II. Third, the CAMRAD II model used a 
low-resolution — just 20 aerodynamic panels span- 
wise and 24 azimuthal stations — representation of 
the rotor disk. This low-fidelity representation 
doubtlessly fails to capture all of the nuances of the 
spanwise loading on the blade. CAMRAD II is ca- 
pable of including many more aerodynamic panels 
in its blade model, and doing so should henceforth 
be considered standard practice. Improving force 
conservation in the coupling process is an area that 
requires future work. 

Since the OVERFLOW 2 model contains the 
most complete representation of the rotor aero- 
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Figure 5 - Thrust Trim for Speed Sweep. 

dynamics, only its data will be presented in the 
remainder of this paper. The reader should bear in 
mind that all results presented here are the result 
of a trim solution with a slightly higher (~ 2.5 To) 
than prescribed thrust. 

Recall that the CFD calculations for the speed 
sweep were performed by matching the wall cor- 
rected shaft angles from the test data rather than 
trimming propulsive force. Figure 6a shows the 
difference between corrected and uncorrected shaft 
angle as a function of advance ratio. Figure 6b 
depicts the corresponding propulsive force at each 
speed. Propulsive force is underpredicted at low 
speed and overpredicted at high speed. A portion 
of this prediction error at high speed is due to in- 
accurate trim. The large shaft angles at high speed 
divert a substantial portion of the excess thrust due 
to trim error into propulsive force. 

Figure 7 compares measured and predicted 
power. Figure 7a shows total power which is 
underpredicted by approximately 2%. Figure 7b 
is formed by removing the parasite power from 
Fig. 7a leaving the sum of profile and induced 
power, the two of which are inseparable in CFD 
simulation. From Fig. 7b it is clear that the sim- 
ulation is having difficulty predicting profile and 
induced power accurately at the high speed points. 
This is masked in Fig. 7a by a fortuitous overpredic- 
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Figure 6 - Prescribed Shaft Angle and Resultant Propulsive Force for Speed Sweep. 
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Figure 7 - Speed Sweep Power. 
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tion of propulsive force, and hence parasite power, 
at these points. A more accurate trim would result 
in a small thrust reduction, causing a complemen- 
tary reduction in induced power, further degrading 
the correlation throughout Fig. 7 slightly. 

Also shown in Figs. 6 and 7 are predictions us- 
ing the simulated tunnel walls for a single condi- 
tion (fi = 0.15). The C x data from this model has 
been corrected using the same Prandtl-Glauert wall 
correction used for the test data and is included in 
Fig. 6b. Predictions with the tunnel walls underpre- 
dict propulsive force and overpredict power at this 
advance ratio to a greater extent than free air sim- 
ulations. One possibility is that the inviscid walls 
of the simulation necessitate a different wall correc- 
tion than the test measured data. 

Sectional Comparisons 

In Fig. 8 are presented the sectional loads — 
M 2 c n , M 2 c m , M 2 c c — at several radial stations for the 
lowest speed in the speed sweep, ^ = 0.15. All data 
are presented with means included. This is at a 
slightly higher thrust than, but otherwise is similar 
to, the oft-studied C85 13 flight condition from the 
UFI-60 flight test [3,14]. 

Several predicted curves are plotted against the 
test data. The first, in green, was computed us- 
ing the 26.2 million point baseline grid. The sec- 
ond was computed with the 50.7 million point fine 
grid, in red. Although there are instances where 
the fine grid improves on the baseline solution (e.g. 
M 2 c c on the advancing side at 7 r= 0.865), there are 
also instances where the solution on the baseline 
grid is closer to the measured data (e.g. M 2 c m in 
the second quadrant at r /R=o.92). Otherwise, these 
two curves are very much in agreement with each 
other. Because the baseline grid performed so well 
for this low-speed flight condition — where accurate 
capture of the wake is most important — it was de- 
cided to use the baseline grid for the remainder of 
this study. 

In general, the trends in normal force are well 
matched. The normal force pulses due to blade vor- 
tex interaction (BVI) begin to appear at r /R = 0.775 
on both the advancing side (near 1/7 = 70°) and re- 
treating side (near 1/7 = 280°). These BVI events are 
visible in the CFD wake visualization of Fig. 9. The 
magnitudes of the normal force pulses in Fig. 8 are 
accurately captured by the simulation save for a 
small underprediction on the advancing side BVI 
at '/r = 0.865. The largest missed predictions are 


mean shifts at Vr = o. 4 and 0.865. T n addition, a 
small surge in normal force from 1/7=120° to 180° 
is present in the test data for all of the outboard 
stations, but the analysis prefers a nearly continu- 
ous decrease in normal force throughout the second 
quadrant. 

The pitching moment is generally mean shifted 
with a tendency toward overprediction. It should 
be noted that the pitching moment measured dur- 
ing the test is very sensitive to transducer perfor- 
mance at the trailing edge, so a mean shift is likely 
to include significant measurement error. One com- 
mon exception to the trend of overpredicted mean 
occurs at r /R= 0.775. A possible cause is the trim 
tab which extends from 7)1=0.73 to 0.86. The trim 
tab on the tested blades was deflected upward by 
approximately 2° which increases the local pitch- 
ing moment significantly compared to adjacent 
stations. Although the trim tab is represented in 
the CFD grid, it is not deflected, and the predicted 
pitching moment mean changes little compared to 
7 'r=o. 4 or 0.865. The moment pulses due to BVI 
are predicted accurately in phase but the peak-to- 
peak amplitude is generally smaller compared to 
measured data. 

Chord force displays a significant mean shift 
with a tendency toward underprediction. There is 
an additional offset that further reduces the chord 
force only on the advancing side. The sign conven- 
tion for chord force is positive toward the leading 
edge so an underprediction represents an excess 
of drag. One contributor to this disparity between 
measurement and prediction is viscous force. The 
test data is derived from integrated pressure mea- 
surements with no correction for viscous effects, 
while the predicted data is the combination of inte- 
grated pressure data and viscous drag at the airfoil 
surface. Predicted airloads computed from pres- 
sure data alone are also presented in Fig. 8. Note 
that the viscous forces make significant contribu- 
tions solely to chord force, and that they make the 
largest contribution on the advancing side near the 
tip. Flowever, the viscous forces account for at best 
20-30% of the difference between measurement 
and prediction. Despite the mean shift, the trend 
in chord force is well represented. Drag excursions 
coincident with the BVI events are captured by 
analysis, particularly on the retreating side where 
the measured drag pulse is closely matched in 
magnitude, phase, and shape. 

Figure 10 contains the same measured and 
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Figure 8 - Sectional Airloads at fi = 0.15, Cj/a= 0.093. 
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Figure 9 - Wake Visualization From CFD. Isosur- 
face of Q-Criterion With Blades Colored 
By Pressure. 


free-air predicted data from Fig. 8 but compares 
it to data predicted using the tunnel wall model. 
The largest differences between the free-air model 
and the tunnel model appear on the advancing 
side, particularly near the 1/1 = 90° BVI. This event 
appears slightly attenuated compared to the free- 
air model in all three coefficients and at all of the 
stations where it is to be observed. 

The low speed case is when the tunnel walls ex- 
ert their maximal effect on rotor performance. Yet 
the free-air simulation compares well to the cor- 
rected test and tunnel simulation data for both in- 
tegrated and sectional loads. This suggests that the 
simple Prandtl-Glauert correction for angle of at- 
tack is reasonable. 

Figure 11 presents the sectional airloads at the 
highest advance ratio in the speed sweep, }i = 0.4. 
At the two outboard stations the analysis predicts 
a dynamic stall cycle in the fourth quadrant. This 
is identifiable as dips in all three coefficients just 
after 1/1 = 270° . Although there is some evidence of 
stall in the test data (particularly in chord force), 
the magnitude of the stall event was predicted to 
be much larger than was measured. The normal 
force trend is well represented, except for a ten- 
dency to overpredict the negative lift on the ad- 
vancing side. Pitching moment predictions show 
larger differences. From r /R= 0.775 outboard, a mo- 


ment pulse due to supercritical flow appears on the 
advancing side at 1/1 = 90°, but its magnitude is di- 
minished in prediction. For both of the outboard 
stations, the analysis underpredicts the magnitude 
of the moment on the advancing side. Chord force 
results are well correlated at the inboard station, 
r /R=o.4. At r /R= 0.775 the chord force trend on the 
advancing side is completely missed, but this may 
be an effect of the difference in trim tab deflec- 
tion. At the two outboard stations, chord force re- 
mains underpredicted on the advancing side, but 
the shape of the predicted chord force curve fol- 
lows that of the test for the first two quadrants, 
even replicating the two small peaks in the second 
quadrant. Differences in the fourth quadrant are 
consistent with excessive stall in the simulation. 

Figure 12 presents a side-by-side comparison 
of test and predicted normal force coefficient at 
r /R= 0.92 for every speed in the sweep for which a 
simulation was run. The test data exhibits vortex 
induced loading at if) = 70 0 and 280° for the lowest 
advance ratios. At higher advance ratios, the in- 
creasingly negative loading in the second quadrant 
balances increased normal force over the rest of the 
rotor disk. The second quadrant minimum starts at 
approximately 1/1 = 135 0 and draws toward 1/7 = 90° 
as advance ratio increases. 

Figure 13 makes the same comparison as Fig. 12 
but for moment coefficient. In the test data, vortex 
effects are again evident at the lowest advance ra- 
tios near 1/1 = 70° and 280°. A large moment pulse 
consistent with supercritical flow develops near 
1/1 = 80° for the highest advance ratios. The moment 
minimum in the second quadrant migrates from 
1/1 = 135° at }i = 0.15, toward 1/1 = 90° at }i = 0.3, and 
then back to 135° by fi = 0.4. Lastly, the test data 
appears to exhibit a weak stall event near 1/1 = 315° 
at ji = 0.4. 

Observing the evolution of the advancing side 
normal force trough in Fig. 12 reveals that both the 
test and analysis display the same monotonic re- 
duction in normal force with speed. The change in 
minimum normal force from one speed to the next 
appears slightly larger for the prediction than in 
test. The simulation also predicts significant lift loss 
from stall in the fourth quadrant at the two highest 
speeds. No corresponding feature is present in the 
measured data. Despite the differences in moment 
coefficient magnitude between test and prediction 
in Fig. 13, the trend between speeds is very sim- 
ilar. The advancing side shock grows in strength 
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Figure 10 - Sectional Airloads at }i = 0.15, Cj/ a- 0.093. Free Air Vs. Tunnel Model. 
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Figure 11 - Sectional Airloads at }i- 0.4, Cj/ u= 0.093. 
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Figure 12 - Normal Force Coefficient at r /R=o.92 for }i = 0.15, 0.2, 0.3, 0.37, 0.4. Cj/cr-o.ogj. 
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Figure 13 - Moment Coefficient at r /R=o.92 for pi = 0.15, 0.2, 0.3, 0.37, 0.4. Cj/a-o.og^. 
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and migrates toward ip = go°. The moment trough 
behind the shock also evolves very similarly across 
the speed range. Differences include an expansion 
of the trough into the third quadrant for the pre- 
dicted data, the prediction of a stall cycle at }i = 0.37 
when none is present in measured data, and an 
enhancement of the predicted stall cycle at /( = 0.4 
compared to measurement. 

Thrust Sweep 

Integrated Results 

Presented in Fig. 14 are the predicted perfor- 
mance results for the selected thrust sweep. From 
this figure it is clear that the analysis performs well 
at predicting the maximum rotor thrust and the re- 
lationship between power and thrust throughout 
the sweep. Recall that the analysis was trimmed 
to match thrust at the Ct/(7=o.o8 point but that 
in general CFD integrated thrust is slightly higher 
than specified. This leads to a small thrust surplus 
throughout the linear region since all collective an- 
gles are set relative to the baseline condition. At 
high collectives the excess thrust vanishes because 
the rotor is stalled and appears incapable of pro- 
ducing more thrust. The analysis also exhibits re- 
duced sensitivity to collective change evidenced by 
an approximately 7% reduction of slope in the lin- 
ear region. 

Power is generally overpredicted, with differ- 
ences up to 11% at the highest thrust conditions 
(Fig. 14b). The fact that the predicted and measured 
points lie on very similar curves, however, suggests 
that if thrust was rigorously matched, power may 
be better predicted. 

Sectional Comparisons 

Sectional airloads for the minimum collective 
setting are presented in Fig. 15. At this collective, 
the measured thrust was Qj/ a =0.02. The normal 
force bucket in the second quadrant is narrower 
and shallower in prediction than in test. Normal 
force is also slightly overpredicted for the majority 
of the retreating side, except at the inboard station. 
Moment coefficient mean is well predicted at this 
condition — except at r /R=0-4 — but there exist sig- 
nificant differences in trends between analysis and 
test. For example, there is a shock-induced mo- 
ment pulse that appears in the second quadrant 
at r /R=0-775, migrates toward the first quadrant. 


and ends in an interesting double pulse at r /R=o.92. 
Although the phase of these pulses is well corre- 
lated, the peak-to-peak magnitude of the predicted 
pulses is smaller. Moment recovery following the 
shocks also begins too early and progresses slowly 
compared to measurement for r /R= 0.865 an ^ 0.92. 
Chord force continues to be mean shifted outboard 
of r /R= 0.775, a n d the trend is completely mispre- 
dicted at r /R= 0.775. 

The sectional airloads presented in Fig. 16 
are for the maximum collective condition where 
Cx/u=o.i255. Comparisons at this condition must 
be tempered with the knowledge that stall is by 
nature a chaotic phenomena and there exists con- 
siderable variation in loading from one revolution 
to another, even in simulation. Test data in Fig. 16 
are average loads sampled from 128 revolutions, 
but computing an ensemble average from CFD is 
impractical. The predicted data provided in Fig. 16 
are instead single-revolution snapshots. 

The flow at the inboard station remains rela- 
tively benign and correlation is good. The most 
significant deviations are a familiar mean shift in 
pitching moment and a blunting of the chord force 
peak in the second quadrant. Outboard, the load- 
ing is dominated by two severe stall cycles clearly 
indicated near the forth quadrant. The analysis ap- 
pears to phase shift the stall cycles slightly, pre- 
dicting both events perhaps 15 0 earlier than mea- 
sured in test. Despite the phase shift, the magni- 
tudes of the stall events are well predicted. One 
glaring difference at this deeply stalled condition 
is a large first-quadrant excursion in all three coef- 
ficients for the prediction. The similarity between 
the frequency of this event and the torsional mode 
of the blade suggests that it is a torsional response 
to the fourth quadrant stall cycles. 

Figure 17 compares sectional normal force at 
Vr— 0.92 for all collective settings simulated. The 
experiment observed an interesting evolution of 
normal force in the second quadrant. The normal 
force increased with collective until Cj / cr=o.i2 af- 
ter which it began to decrease, forming a narrow 
trough at approximately ip = 120°. The normal force 
reduction beginning at Cj / c = 0.1 2 is necessary to 
maintain roll balance across the rotor. Also visi- 
ble in the test data is the development of two stall 
cycles, one at ip = 290 0 and the second at 340°. 

Figure 18 presents moment coefficient at r /R=o.92 
for all collective settings simulated. The most 
prominent features in the test data are the two 
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Figure 14 - Integrated Performance Results for Collective Sweep. n = 0.3, M tip = 0.625, a s = o°. 


moment stall events that develop in the fourth 
quadrant. Both stall cycles appear abruptly at 
Ct/ ( 7=0.12. The first stall cycle quickly achieves 
a stable minimum moment which slowly drifts 
toward = 270°. This is in contrast to the second 
stall cycle where the moment decreases rapidly and 
monotonically but remains fixed in azimuth. 

Changes in normal force from one collective 
setting to the next are similar between test and 
prediction, with only a few notable differences. 
Foremost are the advancing side oscillations in the 
predicted data. This oscillation begins abruptly 
at Ct / <J— 0.120 and increases in magnitude along 
with collective. Examining the predicted moment 
coefficient data in Fig. 18 reveals that there is a sin- 
gle weak stall event at Cj/cr=o.io which quickly 
strengthens and spawns a second stall event by 
Ct /c 7 ”— 0.120. This suggests that the advancing 
side oscillation present in the predicted data is a 
response to either the strengthening of the first 
stall event, the formation of the second stall event, 
or both. Figure 19 plots the local pitch angle at 
t /r= 0.92 for three of the highest collective settings. 
At Cj/a-o.io the rotor has a single stall cycle and 
very little torsional oscillation is visible in Fig. 19. 
As thrust increases, the fourth quadrant stall cy- 
cles strengthen concurrently with a 5 / rev torsional 
oscillation in Fig. 19. This 5/ rev oscillation grows 


in magnitude and ultimately persists through all 
four quadrants at Ct/ £ 7=0.1 255. Another differ- 
ence is the second quadrant negative loading in 
Fig. 17 which is distorted in the predicted data. 
The trough is diminished in magnitude and shifted 
toward the end of the second quadrant. This differ- 
ence appears to be a consequence of the torsional 
oscillation persisting into the second quadrant. 

The effects of the torsional oscillation are also 
readily visible in the predicted data for Fig. 18. 
The oscillation dominates the predicted data in 
the first quadrant for the high thrust cases. The 
oscillation may also be responsible for the missed 
prediction of the moment minimum appearing in 
the test data at ip = 135 0 . Examine the stall cycles 
in the fourth quadrant. The predicted responses 
are more chaotic and appear to have higher fre- 
quency content than the measured data. The first 
cycle marches toward 270° in an orderly fashion 
according to test data, but the analysis predicts a 
wider azimuthal range for the minima and high- 
frequency oscillations make it difficult to discern 
any orderly progression. The final location of the 
first stall cycle — at Cj/cr=o.i 255 — is very similar 
between test and prediction. Comparisons between 
test and prediction for the second stall cycle in- 
clude many of the same observations. The high 
frequency content of the predicted data make it 
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Figure 15 - Sectional Airloads, Cj/cr=o.02. 
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Figure 16 - Sectional Airloads, Ct/V=o.1255. 
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Figure 17 - Normal Force Coefficient at r /R=o.92 for Cj/cr=o.02, 0.04, 0.08, o.io, 0.120, 0.123, 0.1247, 0.1253 
0.1255. (Measured values, predicted Cj/j is slightly higher.) 
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Figure 18 - Moment Coefficient at r /R=o.92 for Ct/c^o. 02, 0.04, 0.08, 0.10, 0.120, 0.123, 0.1247, 0.1253 
0.1255. (Measured values, predicted Cj/ J is slightly higher.) 
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Figure 19 - Local Pitch Angle at r /R=o.92, 2 /rev and 
Above. 

difficult to accurately locate the minimum, but the 
majority seem clustered near 320°. The test shows 
a slightly later location for the second stall cycle, 
around 330°. 

Conclusions 

This paper has described the initial effort to 
correlate prediction with wind tunnel measure- 
ments from the recently completed UH-60A Air- 
loads wind tunnel test. Differences identified in 
the present work represent opportunities to im- 
prove prediction methodologies or identify defects 
in measured data. 

The simulation methodology involved a com- 
prehensive rotorcraft simulation package, CAM- 
RAD II, augmented with the first-principles aero- 
dynamics of a CFD simulation provided by OVER- 
FLOW 2. 

Two parametric sweeps from the wind tunnel 
test were modeled: a speed sweep at constant lift 
and a thrust sweep at constant speed. Simulation of 
both sweeps were complicated by a disagreement 
between the two simulation codes on the integrated 
trim targets. The largest effect of this disagreement 
is a discrepancy of 2.5% in the trimmed thrust be- 
tween the two codes, with the CFD being the higher 
of the two. Because they are a more complete rep- 


resentation of the aerodynamics, only loads from 
CFD have been used for the present work. Further 
work is required to improve the coupling mecha- 
nism and CAMRAD II model to minimize discrep- 
ancies in loading between the two codes. 

For the speed sweep, the analysis displayed a 
tendency to underpredict propulsive force at low 
speed and overpredict it at high speed. This may 
be partially attributable to inaccurate trim. The 
analysis underpredicted total power by approxi- 
mately 2%. Although total power appears to be 
consistently predicted across the speed range, at 
high speed it is just a fortunate result of com- 
pensating errors — overpredicted parasite power on 
one side with underpredicted profile and induced 
power on the other. When applied to the thrust 
sweep, the analysis was very accurate at predicting 
maximum thrust and at capturing the relationship 
between thrust and power. Before the rotor stalled, 
there was a small increase in thrust attributable to 
inaccurate trim in the simulation at the baseline 
condition. The predicted thrust also exhibited a 
reduced sensitivity to collective pitch resulting in a 
shallowing of the slope in the thrust vs. collective 
curve. Power was consistently overpredicted — as 
much as 11% at the peak of stall. However, the vari- 
ation of power with thrust is very similar between 
simulation and test. This suggests that power may 
be better predicted if the simulation trim were more 
accurate. 

Detailed sectional loading comparisons reveal a 
number of differences between measurement and 
prediction. In general, correlation of normal force 
was good, with the prediction only occasionally de- 
viating from test measurement. Comparisons of 
pitching moment exhibited more variation. Pitch- 
ing moment predictions were frequently offset from 
measurement. It was noted that these comparisons 
are hampered by the difficulty of accurately mea- 
suring the mean pitching moment in test. Notewor- 
thy flow features such as blade-vortex interaction, 
supercritical flow, and stall are accurately located 
in azimuth but the predicted magnitude of these 
events is often diminished. Chord force predictions 
bore the least resemblance to test data. These com- 
parisons were often characterized by mean shifts 
and a large overprediction of drag on the advanc- 
ing side of the rotor. As much as 20-30% of this 
difference is attributable to viscous forces included 
in predictions but absent from test data. Correla- 
tion of all three sectional coefficients at r / R = 0.775 


19 


may be improved by correctly modeling trim tab 
deflection. 

Other observed trends included a tendency for 
the analysis to overpredict the severity of stall, fre- 
quently predicting it at lower advance ratio and fur- 
ther inboard on the blade than indicated by mea- 
surement. Individual stall events are not as well 
defined in prediction as in test and their azimuthal 
locations do not progress smoothly as thrust in- 
creases. Also, for flows exhibiting a shock at or 
near ip = go°, the moment trough that formed in the 
second quadrant just after the shock generally ex- 
tended into the third quadrant according to simu- 
lation, but the corresponding feature was confined 
to the second quadrant in the test. 

One final observation was made only for the 
deeply stalled cases of the thrust sweep. The analy- 
sis predicted a large oscillation in all three sectional 
coefficients beginning in the first quadrant and per- 
sisting into the second. No such excursion was 
present in the measured data. This has been shown 
to be a torsional response of the blade model to the 
formation of stall cycles in the fourth quadrant. 

Modeling wind tunnel walls in the CFD simu- 
lation made only small differences in the sectional 
loads at ^ = 0.15. The largest differences were on the 
advancing side near the 1/7 = 90° BVI. Larger differ- 
ences were observed in the integrated performance 
parameters, with the free-air model outperforming 
the tunnel model. Further work is required to de- 
termine if this is a trend or merely an isolated mis- 
prediction. 

The goal of this work was to provide an ini- 
tial correlation with the newest data from the U.S. 
Army/NASA UFI-60 test program. This dataset is 
of tremendous value to the rotorcraft community 
for the breadth of flight conditions sampled and 
the variety of detailed measurements made at each. 
No significant anomalies were identified in the test 
data studied, but the correlation effort revealed a 
number of areas for improvement in the simulation 
technique. Development is required to improve 
the coupling mechanics between CAMRAD II and 
OVERFLOW 2. Better force transfer between the 
two is mandatory for improving trim accuracy and 
accurately predicting integrated loads. Accurate 
power prediction remains a challenge for analysis. 
More detailed investigation is necessary to identify 
the cause of the stall-related torsional oscillation 
in the blade model and to improve the capture of 
sectional loading details such as dynamic stall and 


the effects of supercritical flow. 

References 

[1] Norman, T.R., Shinoda, P., Peterson, R.L., and 
Datta, A., "Full-Scale Wind Tunnel Test of 
the UFI-60A Airloads Rotor", American Fleli- 
copter Society 67 th Annual Forum, May 2011. 

[2] Lorber, P.F., "Aerodynamic Results of a 
Pressure-Instrumented Model Rotor Test at the 
DNW", Journal of the American Helicopter Soci- 
ety, Vol. 36, No. 4, October 1991. 

[3] Kufeld, R.M., Balough, D.L., Cross, J.L., et al., 
"Flight Testing of the UFI-60A Airloads Air- 
craft", American Flelicopter Society 50 th An- 
nual Forum, May 1994. 

[4] Johnson, W., "Technology Drivers in the De- 
velopment of CAMRAD II", American Fleli- 
copter Society Aeromechanics Specialist Meet- 
ing, San Francisco, CA, January 19-21, 1994. 

[5] Yeo, FI., Bousman, W.G., and Johnson, W., 
"Performance Analysis of a Utility Flelicopter 
with Standard and Advanced Rotor", Journal 
of the American Helicopter Society, Vol. 49, No. 3, 
July 2004, pp. 250-270. 

[6] Shinoda, P.M., Yeo, H., and Norman, T.R., "Ro- 
tor Performance of a UFI-60 Rotor System in 
the NASA Ames 80- by 120-Foot Wind Tun- 
nel", Journal of the American Helicopter Society, 
Vol. 49, No. 4, October 2004. 

[7] Yeo, FI., and Johnson, W., "Assessment of 
Comprehensive Analysis Calculation of Air- 
loads on Flelicopter Rotors", Journal of Aircraft, 
Vol. 42, No. 5, September-October 2005. 

[8] Yeo, FI., and Johnson, W., "Prediction of Rotor 
Structural Loads with Comprehensive Analy- 
sis", Journal of the American Helicopter Society, 
Vol. 53, No. 2, April 2008. 

[9] Yeo, FI., Romander, E., and Norman, T., "In- 
vestigation of Rotor Performance and Loads of 
a UFI-60A Individual Blade Control System", 
American Flelicopter Society 66 th Annual Fo- 
rum, Phoenix, AZ, May 11-13, 2010. 

[10] Yeo, FI., Bousman, W.G., and Johnson, W., 
"Performance Analysis of a Utility Flelicopter 


20 


with Standard and Advanced Rotor", Journal 
of the American Helicopter Society, Vol. 49, No. 3, 
July 2004, pp. 250-270. 

[11] Buning, P.G., Gomez, R.J., and Scallion, 
W.I., "CFD Approaches for Simulation of 
Wing-Body Stage Separation", AIAA-2004- 
4838, AIAA 22 nd Applied Aerodynamics Con- 
ference, Providence, RI, August 16-19,2004. 

[12] Potsdam, M., Strawn, R.C., and Meakin, 
R., "Dynamic Rotorcraft Applications Using 
Overset Grids", 31 st European Rotorcraft Fo- 
rum, Florence, Italy, September 13-15, 2005. 

[13] Tung, C., Caradonna, F.X., and Johnson, W., 
"The Prediction of Transonic Flows on an Ad- 
vancing Rotor", American Helicopter Society 
40 th Annual Forum, Arlington, VA, May 16- 
18, 1984. 

[14] Potsdam, M., Yeo, H., and Johnson, W., "Ro- 
tor Airloads Prediction Using Loose Aerody- 
namic/Structural Coupling", Journal of Air- 
craft, Vol. 43, No. 3, May-June 2006. 

[15] Nygaard, T., Saberi, H., Ormiston, R.A., 
Strawn, R.C., and Potsdam, M., "CFD and 
CSD Coupling Algorithms and Fluid Struc- 
ture Interface for Rotorcraft Aeromechanics 
in Steady and Transient Flight Conditions", 
American Helicopter Society 62 nd Annual Fo- 
rum, Phoenix, AZ, May 9-11, 2006. 

[16] Chang, I-C., et al., "Airloads Prediction of a 
UH-60A Rotor Inside The 40- by 8o-Foot Wind 
Tunnel", American Helicopter Society Special- 
ists' Conference on Aeromechanics, San Fran- 
cisco, CA, January 20-22, 2010. 


21 


